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Abstract 

We present the preliminary result of a search for the decay afK*^. The data, collected 

with the BABAR detector at the Stanford Linear Accelerator Center, represent 465 million BB pairs 
produced in e~^e~ annihilation at the T(4S') energy. The result for the branching fraction is: 

B{B+ ^ a+K*^) X B{at ^ 7r+7r-7r+) = (O.?!";^!";^) x 10-^, 

corresponding to an upper limit at 90% confidence level of 1.6 x 10~^. The first error quoted is 
statistical, the second systematic. 
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1 INTRODUCTION 



Recent searches for decays of B mesons to final states with an axial-vector meson ai or bi and a 
pion or kaon have revealed modes with branching fractions that are rather large among charmless 
decays: (15-35) x 10"^ for B ai{TT,K) P E], and (7-11) xlO~^ for a charged pion and kaon 
in combination with a 6^ ^ ^i^ meson [3l H]. On the other hand the experimental search for 
B^ p"^ set an upper limit of 1.7 x 10~^ at the 90% confidence level for the branching fraction 
|5J, although a branching fraction of 25 x 10~^ was expected [6|. 

The available theoretical estimates of the branching fractions of B'^ mesons to afK*^ come 
from calculations based on naive factorization [7], and on QCD factorization [6]. The latter incor- 
porates light-cone distribution amplitudes evaluated from QCD sum rules, and predicts branching 
fractions in quite good agreement with the measurements for B — > aivr and B — > aiK [H [2]. The 
expected branching fraction for B~^ afK*^ from naive factorization is 0.51 x 10~^ and that from 
QCD factorization is 9-7l35l2^4^ ^ 10~® with a prediction for the longitudinal polarization fraction 
fi = 0.38]'^Q4j. The first theoretical error corresponds to the uncertainties due to the variation 
of Gegenbauer moments, decay constants, quark masses, form factors and a B meson wave func- 
tion parameter. The second theoretical error corresponds to the uncertainties due to the variation 
of penguin annihilation parameters [B]. For the longitudinal polarization fraction, all errors are 
added in quadrature, since the theoretical uncertainty is dominated uncertainties in the size of the 
penguin-annihilation amplitude. This mode is expected to be substantially enhanced by penguin 
annihilation and thus it is important to study this mechanism. In this paper we present the first 
search for the decay B^ a^K*^. 

2 THE BABAR DETECTOR AND DATASET 

The data for this measurement were collected with the BABAR detector [8j at the PEP-II asymmetric 
e~^e~ collider located at the Stanford Linear Accelerator Center. An integrated luminosity of 424 
fb^^, corresponding to (465 it 5) x 10^ BB pairs, was produced in e"^e~ annihilation at the 
T(45) resonance (center-of-mass energy ^/s = 10.58 GeV). Charged particles from the e"'"e~ 
interactions are detected, and their momenta measured by a combination of five layers of double- 
sided silicon microstrip detectors and a 40-layer drift chamber both operating in the 1.5 T magnetic 
field of the BaBar superconducting solenoid. Photons and electrons are identified with a CsI(Tl) 
electromagnetic calorimeter (EMC). Further charged particle identification (PID) is provided by 
the average energy loss (dE/dx) in the tracking devices and by an internally reflecting ring imaging 
Cherenkov detector (DIRC) covering the central region. 

A detailed Monte Carlo program (MC) is used to simulate the B meson production and 
decay sequences, and the detector response [9j. Dedicated signal MC events for the decay 5+ — > 
afK*^ with af — > p^n^ has been produced. For the ai(1260) meson parameters we use a mass 
of 1230MeV/c^ and a width of 400MeV/c^. We account for the uncertainties of these resonance 
parameters in the determination of systematic uncertainties. The — > vr^vr^vr^ decay proceeds 
mainly through the intermediate states (7r7r)p7r and (7r7r)o-7r [10]. No attempt is made to separate 
contributions of the dominant P wave (vrTr)^ from the S wave (vr7r)o- in the channel vrvr. The 
difference in efficiency for the S wave and P wave cases is accounted for as a systematic error. 
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3 ANALYSIS METHOD 



We reconstruct a'^ candidates through the decay sequence af p tt^ and p — > vr+vr". The other 
primary daughter of the B meson is reconstructed as K*^ — > K^tt~. For the p^, the invariant 
mass of the pion pair is required to he between 0.55 and 1.0 GeV/cP. For the ai and K* we accept 
a range that includes sidebands. The oi invariant mass is required to he between 0.9 and 1.8 
GeV/c^, while the K* invariant mass is required to lie between 0.8 and 1.0 GeV/c^. Secondary 
charged pion candidates from oi and K* decays are rejected if classified as protons, kaons, or 
electrons by their DIRC, dE/dx, and EMC FID signatures. We reconstruct the B meson candidate 
by combining the four- momenta of a pair of primary daughter mesons, using a fit that constrains 
all particles to a common vertex. From the kinematics of T(4S') decay we determine the energy- 
substituted mass TTT-Es = — and energy difference AE = Eb — ^a/s, where {Eb,Pb) is the 
B meson four-momentum vector, and all values are expressed in the T(4S') rest frame. We require 
5.25 GeV/c^ < tties < 5.29 GeV/c^ and lAii^l < 100 MeV. To reduce fake meson candidates we 
require a B, ai and K* vertex probability > 0.01. 

We also impose restrictions on the helicity-frame decay angle 9k* of the K* mesons. The helicity 
frame of a meson is defined as the rest frame of the meson with the z axis along the direction of 
boost to that frame from the parent rest frame. For the decay K* Kn, Ok* is the polar angle 
of the daughter kaon, and for ai — > /jvr, is the polar angle of the normal to the ai decay plane. 
We define TLi = cos(^j), where i = {K*,ai). Since many background candidates accumulate near 
\nK* \ = 1, we require -0.98 < Hk* < 0.8. 

Backgrounds arise primarily from random combinations of particles in continuum e~^e~ — > qq 
events {q = u,d, s,c). We reduce these with a requirement on the angle between the thrust axis 
[llj of the B candidate in the T{4S) frame and that of the charged tracks and neutral calorimeter 
clusters in the rest of the event (ROE). The distribution is sharply peaked near |cos0t| = 1 
for jet-like continuum events, and nearly uniform for B meson decays. The requirement, which 
optimizes the expected signal yield relative to its background-dominated statistical uncertainty, is 
I cos I < 0.8. BB background arising from b ^ c transitions is suppressed by removing events 
with D meson candidates, reconstructed in the decays — > K'tt^ and K^tt^tt'^, with an 

invariant mass within ib0.02GeV/c^ of the nominal mass value. 

The number of events which pass the selection is 15802. Besides the signal events, these samples 
contain qq (dominant) and BB with b ^ c combinatorial backgrounds, and a fraction of other 
charmless BB background modes. The average number of candidates found per event in the selected 
data sample is 1.5 (2.0 to 2.4 in signal MC depending on the polarization). We choose the candidate 
that is most likely a signal decay, as judged from the output of a neural network, where we use the 
p meson mass and the vertex fit probabilities of B, ai and K* candidates as input variables. 

We discriminate further against qq background with a Fisher discriminant [12] that combines 
four variables: the polar angle of the B candidate momentum and of the B thrust axis with respect 
to the beam axis in the T{4S) rest frame; and the zeroth and second angular moments -Lo,2 of 
the energy fiow, excluding the B candidate, with respect to the B thrust axis. The moments are 
defined by Lj = J2iPi^ |cos6'j|-' , where 6i is the angle with respect to the B thrust axis of a track 
or neutral cluster i, pi is its momentum. 

We obtain yields and the longitudinal polarization from an extended maximum likelihood 
(ML) fit with the input observables AE, mES; the resonance masses and rriK* and the 
helicity distributions TCk* and Tia^ . 

Since the correlation between the observables in the selected data and in MC signal events is 
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small, we take the probability density function (PDF) for each event to be a product of the PDFs 
for the individual observables. Corrections for the effects of possible correlations are made on the 
basis of MC studies described below. 

We determine the PDFs for the signal and BB background components from fits to MC samples. 
We develop PDF parameterizations for the combinatorial background with fits to the data from 
which the signal region (5.26 GeV/c^ < twes < 5.29 GeV/c^ and |Ai?| < 60 MeV) has been 
excluded. 

The TTT-Es and AE distributions are parametrized as linear combinations of the so-called Crystal- 
Ball function [13] and Gaussian. In case of rriEs for qq and BB background we use the threshold 
function x\/l — exp [—^(1 — 2:^)], with argument x = 2mEs/\/s a-nd shape parameter ^. This 
function is discussed in more detail in Ref. [13] • In case of AE for qq and BB background we use 
a polynomial function. The PDFs for the Fisher discriminat Vj{J-) are parametrized as single or 
double Gaussian. The PDFs for the invariant masses of the ai and K* mesons are constructed 
as linear combinations of relativistic Breit Wigner and polynomial functions. We use a joint PDF 
Vj{Ti.K* ,T~(-ai) for the helicity distributions, the signal component is parametrized as the product of 
the ideal angular distribution in TCk* and TCa^ from Ref. [15], times an empirical acceptance function 
G{'Hk* ,'Hai) while the helicity PDF for the other components is the product of the helicity PDFs 
for Ti-K* and Tia^ ■ The Tii distributions in case of qq and BB background are based on Gaussian and 
polynomial functions. We allow the most important parameters (first coefficient of the polynomial 
function for AE, the invariant masses of the ai and the K*, and the width of the Breit Wigner 
for the invariant masse of the K*) for the determination of the combinatorial background PDFs to 
vary in the fit, along with the yields for the signal and qq background. 

The likelihood function is 

C = expl-J2YAl[J2Y,x (1) 

\ i J i j 
V, (niEslV, ir)V, {AE^)V, (ml^ )V, {m\, )V, {Kj,. , K J , 

where A'^ is the number of events in the sample, and for each component j (signal, qq background, 
6 — > c BB background, or charmless BB background), Yj is the yield of component j and Vj{x^) 
is the probability for variable x of event i to belong to component j. 

We validate the fitting procedure by applying it to ensembles of simulated experiments with 
the qq component drawn from the PDF, and with embedded known numbers of signal and BB 
background events randomly extracted from the fully simulated MC samples. By tuning the number 
of embedded events until the fit reproduces the yields found in the data, we find a positive bias 

Table 1: Summary of results for B^ —>■ afK*^. Signal yield Y, fit bias Yf,, the branching fraction 
of K+vr- B{K*'^ K+TT-), branching fraction B = B{B+ afK*'^) x B{at 7r+7r-7r+), 

significance S (see text) and upper limit UL. The given uncertainties on fit yields are statistical 
only, the uncertainties on the fit bias include the corresponding systematic uncertainties. 

Y Yb BiK*'> K+TT-) e(10-«) S UL (10"^) 

55li? 27 ±14 I 0.7t[J:^l[]:^ 0.9 1.6 



10 



Yf,, to be subtracted from the observed signal yield Y: the corresponding numbers are reported in 
Table [H 

In the fitting procedure we allow the longitudinal polarization to vary, finding the best value 
fi = 1.1 lb 0.2, where the error is statistical only; systematic uncertainties are not evaluated and 
we do not report a measurement on this quantity, since the observed signal is not statistically 
significant. 

We compute the branching fraction by subtracting the fit bias from the measured yield, and 
dividing the result by the number of produced BB pairs and by the product of the selection ef- 
ficiency times the branching ratio for the B{K*^ — > K~^7r~) decay. We assume that the branching 
fractions of the T(4S') to B^B~ and B^B^ are equal, consistent with measurement |T0]. The effi- 
ciency for longitudinally and transversely polarized signal events, obtained from MC signal model, 
is 12.9% and 18.6%, respectively. The results are given in Table [H along with the significance, 
computed as the square root of the difference between the value of — 21n£ (with additive system- 
atic uncertainties included) for zero signal and the value at its minimum. In order to obtain the 
most conservative upper limit, we assume = 1 in estimating the branching fraction. In Figure [U 
we show the projections of data with PDFs overlaid. The data plotted are subsamples enriched in 
signal with the requirement of a minimum value of the ratio of signal to total likelihood, computed 
without the plotted variable. 

4 SYSTEMATIC STUDIES 

Systematic uncertainties on the branching fractions arise from the imperfect knowledge of the 
PDFs, BB backgrounds, fit bias, and efficiency. PDFs uncertainties not already accounted for 
by free parameters in the fit are estimated from varying the signal-PDF parameters within their 
uncertainties. For K* resonance parameters we use the uncertainties from Ref. |10j and for the ai 
resonance parameters from Ref. jl] . The uncertainty from fit bias (Table [T]) includes its statistical 
uncertainty from the simulated experiments, and half of the correction itself, added in quadrature. 
We vary the BB background component yields by 100% for charmless background and by 20% for 
the b ^ c BB background. 

In the systematic uncertainty we account for a possible B'^ — > a2K*^ contribution by parame- 
terizing its PDFs on a dedicated sample of simulated events; for the helicity part of this component 
we use the corresponding joint ideal angular distribution from Ref. [15], as we do for our signal 
component. We conservatively assume B~^ — > a^K*'' branching ratio could be as large as the 
B~^ — > afK*^ and vary the B~^ — > a2K*^ from to 19 events. 

The uncertainty from the polarization is obtained by varying fi within errors found in studies 
where fi was allowed to vary in the fit. Uncertainties in our knowledge of the tracking efficiency in- 
clude 0.3% per track in the B candidate. The uncertainties in the efficiency from the event selection 
are below 0.6%. We determine the systematic uncertainty on the determination of the integrated 
luminosity to be 1.1%. All systematic uncertainties on the branching fraction are summarized in 
Table M 

5 RESULTS 

We obtain as a preliminary result for the product of branching fractions: 

B{B+ ^ atK*°) X B{ai ^ 7r+7r~7r+) = {O.Jtoitol) x 10-^ 
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Table 2: Summary of systematic uncertainties of the determination of the afK*^ branching 

fraction. 



Source of systematic uncertainty 

Additive errors (events) 

b ^ c BB background 6 

Charmless BB background 12 

5+ a^A'*° background 14 

oi meson parametrization 4 

PDF parametrization 3 

Variation on fi 2 

ML Fit Bias 14 

Total additive (events) 26 



Multiplicative errors (%) 

Tracking efficiency 1.2 

Determination of the integrated luminosity 1.1 

MC statistics (signal efficiency) 0.6 
Differences in the selection efficiency for the ai decay 3.3 

Particle identification (PID) 1.4 

Event shape restriction (cos^t) 1-0 

Total multiplicative (%) 4.1 

Total systematic error [-6(10-^)] ±0.7 



corresponding to an upper limit of 1.6 x 10 . 

Assuming ^(af (1260) vr+vr-Tr^) is equal to ^(af (1260) t:^t:^-k^), and that S(af (1260) 
37r) is equal to 100%, we obtain: 

B{B+ - a+K*0) = (1.5tJ:01|:l) X 10"^, 

corresponding to an upper limit of 3.3 x 10^®. The first error quoted is statistical and the second 
systematic. Since the signal significance is 0.9 standard deviations, we quote a 90% confidence level 
upper limit. 

The upper limit from this measurement is, on the one hand, in agreement with the prediction 
from naive factorization [7j, and on the other hand, significantly lower than the QCD factorization 
estimation [6], though not inconsistent with it. 
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Figure 1: Distributions for signal-enhanced subsets (see text) of the data projected onto the fit 
observables for the decay — > afK*^; (a) m-ES) (b) ^E, (c) J^, (d) m{piT) for the ai candidate, 
(e) m{K'iT) for the K* candidate, (f) Tix* and (g) Tiai- The solid hues represent the results of 
the fit, and the dot-dashed and dashed lines the signal and background contributions, respectively. 
These plots are made with cuts on the ratio of signal to total likelihood. With respect to the 
nominal fit 19% to 46% (depending on the variable) of signal events remain. 
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